Search for: All records

Abstract Land surface temperature (LST) is crucial for understanding earth system processes. We expanded the Advanced Baseline Imager Live Imaging of Vegetated Ecosystems (ALIVE) framework to estimate LST in near‐real‐time for both cloudy and clear sky conditions at a five‐minute resolution. We compared two machine learning (ML) models, Long Short‐Term Memory (LSTM) networks and Gradient Boosting Regressor (GBR), using top‐of‐atmosphere observations from the Advanced Baseline Imager (ABI) on the GOES‐16 satellite against observations from hundreds of observation sites for a five‐year period. Long Short‐Term Memory outperformed GBR, especially at coarser resolutions and under challenging conditions, with a clear sky R²of 0.96 (RMSE 2.31K) and a cloudy sky R²of 0.83 (RMSE 4.10K) across CONUS, based on 10‐repeat Leave‐One‐Out Cross‐Validation (LOOCV). GBR maintained high accuracy and ran 5.3 times faster, with only a 0.01–0.02 R²drop. Feature importance revealed infrared bands were key in both models, with LSTM adapting dynamically to atmospheric changes, while GBR utilized more time information in cloudy conditions. A comparative analysis against the physically based ABI_LSTproduct showed strong agreement in winter, particularly under clear sky conditions, while also highlighting the challenges of summer LST estimation due to increased thermal variability. This study underscores the strengths and limitations of data‐driven models for LST estimation and suggests potential pathways for integrating ML models to enhance the accuracy and coverage of LST products. 

Abstract Top‐down entrainment shapes the vertical gradients of sensible heat, latent heat, and CO₂fluxes, influencing the interpretation of eddy covariance (EC) measurements in the unstable atmospheric surface layer (ASL). Using large eddy simulations for convective boundary layer flows, we demonstrate that decreased temperature gradients across the entrainment zone increase entrainment fluxes by enhancing the entrainment velocity, amplifying the asymmetry between top‐down and bottom‐up flux contributions. These changes alter scalar flux profiles, causing flux divergence or convergence and leading to the breakdown of the constant flux layer assumption (CFLA) in the ASL. As a result, EC‐measured fluxes either underestimate or overestimate “true” surface fluxes during divergence or convergence phases, contributing to energy balance non‐closure. The varying degrees of the CFLA breakdown are a fundamental cause for the non‐closure issue. These findings highlight the underappreciated role of entrainment in interpreting EC fluxes, addressing non‐closure, and understanding site‐to‐site variability in flux measurements. 

We investigate how effective surface length scales (Le f f ) and atmospheric boundary layer stability modulate surface-induced secondary circulations over a realistic heterogeneous sur- face. The evolution of the circulations and their impact on surface-atmosphere fluxes are studied using coupled large eddy simulations of the CHEESEHEAD19 field campaign. The heterogeneity-induced circulations were diagnosed using time and ensemble averaging of the atmospheric fields. Simulations were performed for summer (August) and autumn (Septem- ber) Intensive Observation Periods of the field campaign, characterised differently in terms of normalised surface length scales and ABL stability. Quasi-stationary and persistent cir- culations were diagnosed in the daytime ABL that span the entire mixed layer height (zi ). Their variation in time and space are presented. Homogeneous control runs were also per- formed to compare and contrast spatial organisation and validate the time-ensemble averaging operation. In the convective boundary layers simulated during the summer time simulations, wavelengths that scale as the effective surface heterogeneity length scales contribute the most to the heterogeneity-induced transport. Contributions from surface-induced circulations were lower in the simulated near-neutral BL for the autumn simulations. We find that both Le f f /zi and ABL static stability control the relative contribution of surface-induced circulations to the area averaged vertical transport. This scale analysis supports prior work over the study domain on scaling tower measured fluxes by including low frequency contribution. We believe that the conceptual framework presented here can be extended to include the effects of sub-grid land surface heterogeneity in numerical weather prediction and climate models and also to further explore scale-aware scaling methodologies for near surface-atmosphere exchanges. 

Abstract In the last decades the energy-balance-closure problem has been thoroughly investigated from different angles, resulting in approaches to reduce but not completely close the surface energy balance gap. Energy transport through secondary circulations has been identified as a major cause of the remaining energy imbalance, as it is not captured by eddy covariance measurements and can only be measured additionally with great effort. Several models have already been developed to close the energy balance gap that account for factors affecting the magnitude of the energy transport by secondary circulations. However, to our knowledge, there is currently no model that accounts for thermal surface heterogeneity and that can predict the transport of both sensible and latent energy. Using a machine-learning approach, we developed a new model of energy transport by secondary circulations based on a large data set of idealized large-eddy simulations covering a wide range of unstable atmospheric conditions and surface-heterogeneity scales. In this paper, we present the development of the model and show first results of the application on more realistic LES data and field measurements from the CHEESEHEAD19 project to get an impression of the performance of the model and how the application can be implemented on field measurements. A strength of the model is that it can be applied without additional measurements and, thus, can retroactively be applied to other eddy covariance measurements to model energy transport through secondary circulations. Our work provides a promising mechanistic energy balance closure approach to 30-min flux measurements. 

Abstract Dry deposition is the second largest tropospheric ozone (O₃) sink and occurs through stomatal and nonstomatal pathways. Current O₃uptake predictions are limited by the simplistic big‐leaf schemes commonly used in chemical transport models (CTMs) to parameterize deposition. Such schemes fail to reproduce observed O₃fluxes over terrestrial ecosystems, highlighting the need for more realistic treatment of surface‐atmosphere exchange in CTMs. We address this need by linking a resolved canopy model (1D Multi‐Layer Canopy CHemistry and Exchange Model, MLC‐CHEM) to the GEOS‐Chem CTM and use this new framework to simulate O₃fluxes over three north temperate forests. We compare results with in situ measurements from four field studies and with standalone, observationally constrained MLC‐CHEM runs to test current knowledge of O₃deposition and its drivers. We show that GEOS‐Chem overpredicts observed O₃fluxes across all four studies by up to 2×, whereas the resolved‐canopy models capture observed diel profiles of O₃deposition and in‐canopy concentrations to within 10%. Relative humidity and solar irradiance are strong O₃flux drivers over these forests, and uncertainties in those fields provide the largest remaining source of model deposition biases. Flux partitioning analysis shows that: (a) nonstomatal loss accounts for 60% of O₃deposition on average; (b) in‐canopy chemistry makes only a small contribution to total O₃fluxes; and (c) the CTM big‐leaf treatment overestimates O₃‐driven stomatal loss and plant phytotoxicity in these temperate forests by up to 7×. Results motivate the application of fully online vertically explicit canopy schemes in CTMs for improved O₃predictions. 

Abstract Understanding the controlling mechanisms of soil properties on ecosystem productivity is essential for sustaining productivity and increasing resilience under a changing climate. Here we investigate the control of topsoil depth (e.g., A horizons) on long‐term ecosystem productivity. We used nationwide observations (n = 2401) of topsoil depth and multiple scaled datasets of gross primary productivity (GPP) for five ecosystems (cropland, forest, grassland, pasture, shrubland) over 36 years (1986–2021) across the conterminous USA. The relationship between topsoil depth and GPP is primarily associated with water availability, which is particularly significant in arid regions under grassland, shrubland, and cropland (r = .37, .32, .15, respectively,p < .0001). For every 10 cm increase in topsoil depth, the GPP increased by 114 to 128 g C m⁻² year⁻¹in arid regions (r = .33 and .45,p < .0001). Paired comparison of relatively shallow and deep topsoils while holding other variables (climate, vegetation, parent material, soil type) constant showed that the positive control of topsoil depth on GPP occurred primarily in cropland (0.73, confidence interval of 0.57–0.84) and shrubland (0.75, confidence interval of 0.40–0.94). The GPP difference between deep and shallow topsoils was small and not statistically significant. Despite the positive control of topsoil depth on productivity in arid regions, its contribution (coefficients: .09–.33) was similar to that of heat (coefficients: .06–.39) but less than that of water (coefficients: .07–.87). The resilience of ecosystem productivity to climate extremes varied in different ecosystems and climatic regions. Deeper topsoils increased stability and decreased the variability of GPP under climate extremes in most ecosystems, especially in shrubland and grassland. The conservation of topsoil in arid regions and improvements of soil depth representation and moisture‐retention mechanisms are critical for carbon‐sequestration ecosystem services under a changing climate. These findings and relationships should also be included in Earth system models. 

Abstract How convective boundary‐layer (CBL) processes modify fluxes of sensible (SH) and latent (LH) heat and CO₂(F_c) in the atmospheric surface layer (ASL) remains a recalcitrant problem. Here, large eddy simulations for the CBL show that whileSHin the ASL decreases linearly with height regardless of soil moisture conditions,LHandF_cdecrease linearly with height over wet soils but increase with height over dry soils. This varying flux divergence/convergence is regulated by changes in asymmetric flux transport between top‐down and bottom‐up processes. Such flux divergence and convergence indicate that turbulent fluxes measured in the ASL underestimate and overestimate the “true” surface interfacial fluxes, respectively. While the non‐closure of the surface energy balance persists across all soil moisture states, it improves over drier soils due to overestimatedLH. The non‐closure does not imply thatF_cis always underestimated;F_ccan be overestimated over dry soils despite the non‐closure issue. 

Single point eddy covariance measurements of the Earth’s surface energy budget frequently identify an imbalance between available energy and turbulent heat fluxes. While this imbalance lacks a definitive explanation, it is nevertheless a persistent finding from single-site measurements; one with implications for atmospheric and ecosystem models. This has led to a push for intensive field campaigns with temporally and spatially distributed sensors to help identify the causes of energy balance non-closure. Here we present results from the Chequamegon Heterogeneous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors 2019 (CHEESEHEAD19)—an observational experiment designed to investigate how the Earth’s surface energy budget responds to scales of surface spatial heterogeneity over a forest ecosystem in northern Wisconsin. The campaign was conducted from June–October 2019, measuring eddy covariance (EC) surface energy fluxes using an array of 20 towers and a low-flying aircraft. Across the domain, energy balance residuals were found to be highest during the afternoon, coinciding with the period of surface heterogeneity-driven mesoscale motions. The magnitude of the residual varied across different sites in relation to the vegetation characteristics of each site. Both vegetation height and height variability showed positive relationships with the residual magnitude. During the seasonal transition from latent heat-dominated summer to sensible heat-dominated fall the magnitude of the energy balance residual steadily decreased, but the energy balance ratio remained constant at 0.8. This was due to the different components of the energy balance equation shifting proportionally, suggesting a common cause of non-closure across the two seasons. Additionally, we tested the effectiveness of measuring energy balance using spatial EC. Spatial EC, whereby the covariance is calculated based on deviations from spatial means, has been proposed as a potential way to reduce energy balance residuals by incorporating contributions from mesoscale motions better than single-site, temporal EC. Here we tested several variations of spatial EC with the CHEESEHEAD19 dataset but found little to no improvement to energy balance closure, which we attribute in part to the challenging measurement requirements of spatial EC.